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Crystallographic studies of channel proteins have provided insight into the molecular mechanisms of ion 
channels, even though these structures are obtained in the absence of the membrane and some structural 
portions have remained unsolved. Here we report the gating structure of the membrane-embedded KcsA 
potassium channel using atomic force microscopy (AFM). The activation gate of the KcsA channel is located 
on the intracellular side, and the cytoplasmic domain was truncated to clear the view of this location. Once 
opened, the individual subunits in the tetramer were resolved with the pore open at the center. Furthermore, 
AFM was able to capture the previously unsolved bulge helix at the entrance. A molecular dynamics 
simulation revealed that the bulge helices fluctuated dramatically at the open entryway. This dynamic 
behavior was observed as vigorous open-channel noise in the single-channel current recordings. The role of 
the bulge helices in the open gate structure is discussed. 

The potassium channel exists ubiquitously on the cell membrane, and maintains cellular electrical activity 1 . A 
vast array of gene species for the potassium channel has been identified across the whole range of living 
organisms 2 . The crystal structures of several channel types have been solved and the basic architecture was 
elucidated 3 " 7 . The channel structure is functionally divided into the pore domain and the sensory domain. The 
pore domain is responsible for the common channel function, i.e. the selective K + permeation and gating. In 
contrast, the sensory domain, which is located in the cytoplasm or membrane, modulates the gating function of 
the pore domain. In this arrangement, the environmental changes sensed by the various types of the sensory 
domain are transferred to the pore domain so as to initiate channel activity. Thus, the pore domain is the 
canonical structural unit for all of the potassium channels 8 , and elucidating the gating conformation of the pore 
domain is thus crucial for an understanding the channel function. 

The high resolution structure of the pore domain was first revealed for the KcsA potassium channel (Fig. I) 9 , 
and the ion permeation mechanism through the selectivity filter has been reported 10 " 14 . The open conformation of 
the KcsA channel 15 as well as that of several potassium channel types 1617 demonstrated that the crossed bundle 
helices (the M2 helix) at the intracellular end of the pore domain relax, which then allows ions to have access to the 
pore. The full-length structure of the KcsA channel was subsequently resolved 1819 , showing the straight helix 
bundle to be the cytoplasmic domain (CPD). However, these structural data were obtained in the absence of the 
membrane, and several maneuvers, including the introduction of mutations and fixation with antibodies should 
have been taken for crystallization, which would bias the intact structure. 

Here we report the structure of the activation gate in the pore domain of the KcsA channel in the membrane- 
embedded state, and atomic force microscopy (AFM) is shown to be useful for capturing such global conforma- 
tional changes upon gating in the native state. The AFM provides the surface contour, and thus the structural 
parts protruding from the membrane surface can be delineated 20 " 23 . For certain membrane proteins, a bulky 
oligomeric domain, protruding into either the extracellular or intracellular space, has been resolved by AFM 24-26 , 
and the conformational change of the extra-membranous parts is relevant to the sensory process. In contrast, our 
focus is the less bulky, membrane-embedded pore domain with the activation gate open and closed. In the various 



SCIENTIFIC REPORTS | 3 : 1063 | DOI: 1 0.1 038/srep01 063 



1 



Extracellular 



6.3 nm 




Closed Open 

Figure 1 | Crystal structure of the closed (residue 22-124) and open 
(residue 30-117) conformations of the KcsA channel (pdb code: 1K4C 
and 3F5W). The longitudinal length of the channel was calculated from 
the crystal structure as the distance between the C a atoms furthest apart 
along the channel axis plus their atomic size. For the closed conformation 
the unresolved terminal residue (F125) is taken into considerations. The 
surface of the channel is colored for the charge distribution (red for 
negative and blue for positive). 

potassium channels the pore domain shares a common architecture, 
and some degree of homologous gating transition should take place. 
What configuration does the activation gate take? How do ions in the 
bulk intracellular solution gain access to the pore through the open 
gate? 

In the pore domain, including that in the KcsA channel, the activa- 
tion gate is located close to the cytoplasmic end, where the global 
conformational changes manifest 27,28 . To observe the open and 
closed gate of the membrane embedded KcsA potassium channel 
from the cytoplasmic side, the CPD, which is formed by a straight 
helix bundle and hence hinders the view, was truncated. The CPD- 
truncated KcsA channel retains its native functions of pH-dependent 
gating and single- channel activity 29 . Between the pore domain and 
the CPD, there is a structural region known as the bulge helix (res- 
idue 118-135; Fig. I) 18 . The charged residues are clustered there and 
function as a part of the intracellular pH- sensor 29 " 31 . The bulge helix 
is a critical structure that governs the size for the ion accessing 
aperture in the open conformation. In the CPD -truncated channel, 
one half of the bulge helix is retained at the C-terrninal end (1 18-125), 
and we call this region the part of the bulge helix (H bulge ). In the 
crystal structure of the open state, the H bulge has remained unsolved. 

In this study, the CPD -truncated KcsA channel embedded in the 
membrane was the target of structural and functional investigation in 
the single-molecule level. In the AFM image that was obtained, the 
tetrameric open conformation was revealed, and the crystallography- 
cally- unsolved H bu i ge was captured. The structure of the H bu i ge was 
generated by means of a simulation with molecular dynamics (MD), 
by which the dynamic behavior was examined. Finally, the functional 
consequences of the fluctuating structure were studied using single- 
channel current recordings. This integrated single-molecule 
approach using AFM, MD simulation and single- channel current 
recordings afforded a clear picture of the dynamic structure and 
function of the activation gate of the KcsA channel. 

Results 

Images of the KcsA in the membrane. The liposomes reconstituted 
with the CPD-truncated KcsA channel (residue 1-125) were 
extended over the mica surface, and AFM measurements were 
performed. Distinct islands of flat membranes became evident and 
protruding particles were distributed exclusively on the membrane 
surface (Fig. 2). The individual particles observed at pH 7.5 (Fig. 2A) 
and 4.0 (Fig. 2B) were found to be substantially different in size. 
These particles represent the cytoplasmic end of the CPD-truncated 
KcsA channel for the following reasons. In the reconstituted lipo- 
somes, the KcsA channel is oriented in the "right-side out" configur- 
ation (i.e. the cytoplasmic domain faces to the inside vesicle). This 



orientation has been confirmed by functional measurements using 
liposome-patch and the planar lipid bilayer methods exploiting the 
location of the pH -dependent sensor of the KcsA channel on the 
cytoplasmic side 30 ' 32 . Thus, the KcsA channel was in an upside-down 
orientation (i.e., the cytoplasmic side up) in the flat bilayer. 

KcsA shortens its longitudinal length upon opening. The height 
profile indicates a flat bilayer membrane with a thickness of 
approximately 4 nm on which surface protrusions are present 
(Fig. 2 A, B, lower panels). The histograms show that the height 
from the mica surface, or the longitudinal length of the channel, 
were significantly different at the two different pH values (Fig. 2C): 
6.4 ± 0.3 nm at pH 7.5 (A) and 4.9 ± 0.4 nm at pH 4.0 (B). A 
substantial shortening of the open pore was clearly resolved by AFM. 

While this shortening was expected from the crystal structures of 
the pore domain in the closed and open conformations 915 (Fig. 1), the 
degree of the shortening was not, since the structure of several of 
the residues in the cytoplasmic end in the open conformation has 
remained unsolved. Thus, the AFM image was compared to the 
crystal structure in terms of the longitudinal length of the channel. 
The height for the closed conformation obtained from AFM (6.4 ± 
0.3 nm) matches well with that of the crystal structure (6.3 nm; blue 
arrowhead). In the case of the open conformation, there is a substan- 
tial discrepancy between the height of the AFM image (4.9 nm) and 
that of the crystal structure (4.3 nm; red arrowhead) 9 . This discrep- 
ancy is schematically shown as the light blue square (Fig. 2D), and 
AFM captured the structure of the crystallographically-unsolved C- 
terminal end in the open conformation. 

The open gate and the tetrameric structure. The channels were 
observed at a higher magnification (Fig. 3). The individual par- 
ticles are well resolved for the open channel (Fig. 3B), having a 
half- width of 6.0 nm in the height profile. For the closed channel, 
on the other hand, large particles are seen (Fig. 3A), and the coarse- 
grained particles within them are resolved. The height profile shows 
a periodicity with a half-width of 6.0 nm. This size matches the 
open channel, suggesting that the channel is clustered at a neutral 
pH. 

Further magnification images of the open gate at an acidic pH (4.0) 
are shown (Fig. 3C). The individual subunits in the tetrameric struc- 
ture along with the hole in the center are clearly resolved. In contrast, 
the tetrameric structure was not resolved for the closed conformation 
(Fig. 3 A inset). The six images are snapshots of the open gate 
(Fig. 3C), which would evoke the dynamic picture of the gating. In 
the height profile, a shallow dent is observed in the middle of the 
protrusion (Fig. 3B, the arrow in the lower panel), representing the 
open entrance. These images comprise the first observation of the 
open gate structure of the membrane- embedded KcsA channel. 
Averaging the open gate images reveals the distinct subunit assembly 
around the pore (Fig. 3D). From the height profile of the open con- 
formation, the distance between the subunits was measured for the 
peaks in the subunits. The distance between the diagonal subunits 
was 4.4 ± 0.3 nm (n = 10) and that between the adjacent subunits 
was 3.3 ± 0.3 nm (n = 11). These values are clearly not in agreement 
with those of the crystal structure (3.9 nm and 2.8 nm). 

The evident discrepancies in the molecular height and aperture 
size between the AFM images and the crystal structure indicate that 
the crystallographycally unsolved C-terminal residues, involving the 
H bulge , extend towards the cytoplasmic space. In general, vigorously 
fluctuating structural regions remain unsolved in crystallography. 
Therefore, to obtain a dynamic picture of the AFM-captured struc- 
ture, the terminal regions were reproduced using a computer simu- 
lation method. 

Fluctuation of the iibuige detected by MD simulation. To generate 
the most plausible structure and configuration of the H bu i ge at the 
open pore entrance, the relevant regions were grafted to the open 
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Figure 2 | KcsA images embedded in the membrane. (A, B) AFM images of the CPD -truncated KcsA channel in lipid bilayer and the height profiles. The 
KcsA-reconstituted proteoliposomes were added onto the mica surface, and the liposome membrane was extended on the mica substrate, forming the flat 
bilayer. The cytoplasmic part of the KcsA channel protruded from the membrane surface. The pH of the overlaying buffer solution was 7.5 for A and 4.0 
for B. The height profiles along the green solid lines are shown below the images. Images were recorded with the acoustic AC mode of AFM. (C) Height 
distribution pattern of the CPD -truncated KcsA on the mica surface. The height was evaluated as the sum of the protruded height and the membrane 
thickness. The blue bars indicate the histogram for the closed conformation at pH 7.5 and the red ones the open conformation at pH 4.0. The green lines 
indicate the results of Gaussian fitting. The height was 6.4 ± 0.3 nm at pH 7.5 (n = 170) and 4.9 ± 0.4 nm at pH 4.0 (n = 121). The blue and red arrows 
indicate the height values obtained from the crystallographic structures for the closed (6.3 nm) and open (4.3 nm) states. (D) Schematic representation of 
the membrane-embedded channel on the mica surface. The intracellular structure protruding from the membrane surface was measured with AFM. The 
light blue box shown on the cytoplasmic end of the open conformation represents an AFM-retrieved structural part. 



crystal structures, and an MD simulation of the CPD -truncated 
transmembrane domain (residues 22-125) in the membrane- 
embedded condition was performed. 

The structural configuration of the H bulge in the fully- open struc- 
ture was obtained using a novel method exploiting structural 
templates with various degrees of pore opening. In the closed con- 
formation, the crystallographic data resolved nearly all of the amino - 
acid residues of the CPD-truncated channel up to residue 124. On the 
other hand, the available open structures with various degrees of 
opening left the H bulge unsolved 15 . Starting from the closed structure 
(pdb code: 1K4C), the transmembrane region was made to move 
gradually to the templates of successively more widely open confor- 
mations (pdb code: 3F7Y, 3F7V, 3F5W; see the Methods section). 
The partially neutrarized H bulge at acidic pH (see Methods) was 
grafted at the end of the templates, and the connecting portion 
between the H bu i ge an d the transmembrane region was moved farther 



apart in an centrifugal manner. Meanwhile the H bulge protruding 
towards the cytoplasmic space was allowed to move freely. 

Typical structures in the closed and open states obtained by the 
MD simulation are shown in Fig. 4. The average longitudinal length 
of the simulated closed structure was 6.4 ± 0.05 nm, which is in good 
agreement with both the crystal structure and the AFM images 
(Fig. 4A). In contrast, the longitudinal length of the open structure 
was elongated substantially compared with that shown in both the 
AFM images and the crystal structure. The average length of the 
channel in the open state was estimated to be 5.4 ±0.13 nm, which 
is longer than that in the AFM images by 0.5 nm (shown in magenta 
in Fig. 4C,D). Moreover, the average distance of the diagonal sub- 
units is 4.5 nm, which is slightly more separated compared with that 
in the AFM images (Fig. 4D), while no differences were seen in the 
closed conformation (Fig. 4B). These differences in the open struc- 
ture between the AFM images and MD simulation indicate that the 



SCIENTIFIC REPORTS | 3 : 1063 | DOI: 1 0.1 038/srepOl 063 



3 





30 40 
Distance (nm) 



20 30 
Distance (nm) 




Figure 3 | Structural images of the cytoplasmic end of the membrane-embedded KcsA channel. ( A,B) Cytoplasmic views of the KcsA channel at pH 7.5 
(A) and pH 4.0 (B) and their height profiles along the green line. Particles are resolved in an island at neutral pH, while individual particles are identified at 
acidic pH. (Inset) A 3-times enlarged view of the arrowed KcsA channels. (C) A further magnified view of the open KcsA channel at acidic pH. The scale 
bars indicate 3 nm. (D) An averaged AFM image of the open pore at acidic pH. The averaging was performed from 14 images of the open KcsA channels 
from Figure 3B. The scale bars indicate 3 nm. 



MD simulation complemented the very end of the H bulge in the open 
gate. 

The MD simulation showed that the H bulge sticks out into the 
cytoplasmic space in the open state and fluctuates substantially (in 
yellow and magenta), while the main transmembrane regions (cyan) 
maintain an approximate 4-fold symmetry (see the Supplemental 
video, Movie 1). This flexibility is represented as the B-factor, i.e. 
the temperature factor related to the structural fluctuation 33 (Fig. 4E). 
In both the open and closed states, the B-factors of the transmem- 
brane helices obtained by the MD simulation were coincident with 
that of the closed (1K4C) structure, but smaller than the open 
(3F5W) structure. In contrast, the B-factors of the portions sticking 
out into solution (the extracellular loop and both ends) exhibited a 
higher value due to thermal fluctuations. In particular, the extremely 
high B-factor of the H bulge clearly reflects its fluctuation. On the other 
hand, in the closed state, the cytoplasmic ends form a bundle cross- 
ing, so the fluctuations are suppressed (see the Supplemental video, 
Movie 2). 

The four H bu i ge regions exhibited different structural configura- 
tions, and these were interconvertible with each other on a slower 
time scale (Movie 1). In fact, changes in the structural configurations 
of each of the helices were observed several times over the total 
simulation time of 200 ns, and once a helix obstructed the pore, it 
did so for 5 to 10 ns or more. These results suggest that the move- 
ment of each helix has slower transition kinetics, and the MD simu- 
lation picked up some of these events. 

The fluctuating ff bulge chops up the ionic current. To determine the 
relationship of the structural fluctuation of the H bu i ge observed in the 



MD simulation to the functional behavior of the channel, single- 
channel current measurements were performed. Here we 
examined a non-inactivating mutant of the KcsA channel, E71A, 
without a CPD, such that the functional features of the activation 
gate are demonstrated (Fig. 5) 34 . The single- channel currents of the 
E71A mutant channels at pH 4.0 revealed that the activation gate is 
kept open more than 90% of the time, during which period vigorous 
fluctuations in the current are observed (Fig. 5A). Brief and frequent 
intermissions towards the non-conductive state present, and these 
fast events could not be fully resolved under the limited bandwidth of 
the current recordings (a 2 kHz cut-off frequency). The amplitude 
histograms were drawn with the vertical axis in the logarithmic scale 
such that the current noise was clearly evident 35 ' 36 . Compared to the 
Gaussian distribution (the solid curve) the histogram of the open 
state shows a one-sided shoulder with a long lasting tail, which is 
consistent with the presence of unresolved brief closures 35 ' 37 . A sub- 
level (the red arrow head) is evident in both the current traces and the 
amplitude histogram. 

To address the fluctuating current in relation to the fluctuating 
motion of the H bulge , complete deletion of the H bulge might be helpful. 
However, the deletion would lead to a loss of pH sensitivity, and is 
thus inappropriate for functional measurements. Alternatively, the 
addition of the cytoplasmic domain imposes a load on the freely 
fluctuating bulge helices. Thus, the KcsA channel with and without 
the CPD makes it possible to compare the single channel behavior in 
order to elucidate the role of the CPD and the activity of the H bulge . In 
Fig 5B, the single-channel currents of the full-length channel are 
shown. With the CPD, the current fluctuation was attenuated 
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Figure 4 | Dynamics of the CPD-truncated transmembrane structure of 
the KcsA channel according to the MD simulation. (A) The side view of 
the closed structure (the cytoplasmic side up). (B) The bottom view of the 
closed structure. In the MD -simulated structure, the transmembrane 
region is colored cyan, which is similar to the X-ray crystallographic 
structure. The cytoplasmic end, that missed in the X-ray crystallographic 
structure is shown in yellow. (C) The side view of the open structure. 
(D) The bottom view of the open structure. The transmembrane region is 
colored cyan, which is similar to the X-ray crystallographic structure. The 
cytoplasmic end, which is colored in yellow, constitutes the plausible 
structure that the AFM images detected. The structural parts, colored 
magenta, make up the terminal structure that the MD simulation resolved. 
The membrane was not shown in these figures. (E) The B-factor was 
calculated from the MD simulation. The red and blue lines represent the B- 
factor obtained from the MD simulation for the open and closed 
structures, respectively. The B-factors of the crystal structures in open and 
closed states are also shown (the broken lines). The B-factors for residues 
118-125 are significantly large. The B-factor for the extracellular loop was 
substantially higher than the intra-membrane region. (Inset) A subunit of 
the open configuration is shown, assigning the B-factors for the C a atoms 
with the volume of the spheres. 

substantially, although the open probability was unchanged. The 
shoulder of the open channel in the histogram was depressed sub- 
stantially. 

The more frequent transition to the non-conducting state for the 
CPD-truncated channel indicates that the bulge helices freely fluc- 
tuate and transiently obstruct the conductive pathway. This choppy 
activity is attenuated but not eliminated in the presence of the CPD, 
which restricts the free movement of the bulge helices. The sub-level 



activity occurred exclusively for the CPD-truncated channel, and this 
suggests that the freely moving bulge helices incompletely closed the 
conductive pathway. 

Discussion 

In this study, a dynamic picture of the open gate for the KcsA channel 
embedded in the membrane was obtained by means of an integrated 
single-molecule approach. AFM revealed the novel findings that the 
cytoplasmic end of the transmembrane pore domain protruded from 
the membrane surface and the longitudinal length was substantially 
shortened upon opening. The high- resolution AFM images render 
the clear distinction between the closed and open gate stracture. 
In the open conformation, the tetrameric subunits were separately 
resolved with the pore open at the center. While the longitudinal 
length of the closed conformation was in agreement with that of the 
crystal structure, the substantially elongated and widened aperture 
relative to the crystal structure was also elucidated for the open 
conformation. The AFM -captured structural regions included the 
H bulge , and thus the relevant portion was grafted onto the open 
crystal structure in the course of the MD simulation. The MD simu- 
lation depicted the dynamics of the H bulge in the open gate, which 
manifested as open- channel noise in the single- channel current as 
the result of a "chopping up" of the entryway. 

The time span of the H bulge fluctuation in the MD simulation 
(< hundreds of ns) and the single- channel current recordings 
(> hundreds of |is) differs substantially. The end of the H bulge 
observed in the MD simulation moves rapidly, in the time range of 
ns, and the H bulge moves as a rigid-body, with a correspondingly 
slower time scale. The random hits for each helix and maintaining 
the obstructed configuration for a period of several ns led to the 
collective motions of the helices to slam the open pore shut. This 
stochastic event is likely to be both rare and brief, in a time range of 
|is. On the other hand, the fluctuating single- channel current traces 
suggested events in the time range of sub-ms. However, based on the 
theoretical consideration of the open channel noise 35 , the asymmetric 
shoulder of the amplitude histogram confirms that the underlying 
events occur in a time span of jas. Thus, a gap in the time span for the 
two results in the MD simulation and the single- channel current data 
is merging. 

The open-channel noise in the single- channel current was vig- 
orous in the CPD-truncated KcsA channel (Fig. 5). By restoring 
CPD as a full-length channel, the current fluctuation was substan- 
tially attenuated but still remained. Thus, CPD restricts but does 
not completely eliminate the flexible movement of the H bulge . 
What is the functional role of the fluctuation that takes place in 
the bulge helices? The fluctuating behavior is likely an inherent 
property of the highly charged pH sensor 29,31 which bears net 
positive charges at acidic pH. Earlier reports demonstrated that 
the pH- dependence was reduced by deleting the CPD (The Hill 
coefficient of the pH sensitivity was shown to be 0.7 for the CPD- 
truncated channel, rather than 1.4 for the full-length channel 29 ). 
Based on these findings, we hypothesize that there is a functional 
role of the bulge helix. At neutral pH, the bulge helices are stably 
bundled through either electrostatic interaction or hydrogen-bond 
formation. At acidic pH, the net positive charges act as a repulsive 
force in the bulge helices 29 , leading to their fluctuation. In the 
absence of CPD, vigorous fluctuations in each helix may drive 
the opening of the activation gate in each subunit. On the other 
hand, in the full-length channel with the CPD present, and hence 
restricting the free motion of the bulge helices, the individual 
movements of the bulge helices cooperate so as to induce a con- 
certed effect on the activation gate that maximizes pH- sensitivity. 

Efforts are currently underway to gain insight into the continuous 
spectrum of the structural and functional dynamics of channel 
proteins. 
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Figure 5 | Single-channel current at +100 mV without (A) and with (B) 
CPD at pH 4.0. The all-point amplitude histograms are shown with the 
vertical axis in the logarithmic scale, and the peaks are fitted with the 
Gaussian function. The shoulder of the open conductance peak is 
prominent for the CPD-truncated channel, indicating the open-channel 
noise from the highly fluctuating gating. A low conductance sub-level (the 
red arrowhead) was seen for the CPD-truncated channel. 

Methods 

Materials. The KcsA channel was purified as described previously 38 . Cytoplasmic 
domain (CPD) -truncated KcsA was prepared with previously reported procedure 27 . 
Phospholipid and detergent used in this study are L-a-Phosphatidylcholine Type 
IV-S (asolectin; SIGMA) and n-Dodecyl-(3-D-maltoside (DDM; Dojindo). 

Reconstitution of KcsA into liposome. The phospholipid was dissolved in 
chloroform and then the solvent was evaporated under a stream of nitrogen. The 
resulting lipid film was dried in vacuo for at least 6 h. The dry thin film was hydrated 
with a buffer solution (10 mM HEPES/NaOH (pH 7.5), 200 mM KC1), thereby 
providing a liposomal suspension. A DDM solution was added to the suspension to 
give a co-micellar solution (0.06 wt% DDM of the final concentration). The solution 
of the KcsA dissolved in a buffer (0.06 w% DDM in 10 mM HEPES/NaOH (pH 7.5), 
200 mM KC1) was added to this co-micellar solution (in the range of 3/1-5/1 of lipid/ 
protein ratio (w/w)), and incubated for at least 1 h at 4°C After incubation, the 
co-micellar solution was dialyzed against detergent -free buffer for 30 h and stored at 
4°C The proteoliposomal suspension was used for experiments within 1 week. 

Liposome extension on the mica surface and AFM observation. Flat lipid bilayers 
on a mica surface were prepered with vesicle fusion method 39 ' 40 . The proteoliposomal 
suspension was put onto a mica substrate after the addition of MgCl 2 solution. After 
~15 min of incubation, the solution was gently replaced with buffer at a pH of 4.0 
(10 mM succinic acid pH 4.0, 300 mM KC1) or pH 7.5 (10 mM HEPES pH 7.5, 
300 mM KC1). Images were taken with a Picoplus 5500 (Molecular Imaging) in the 
acoustic AC mode at room temperature. Cantilevers used were BL-AC40TS (tip 
radius: <10 nm, spring constant: —0.1 N/m, OLYMPUS). All images were recorded 
in buffer solution in liquid cells. To evaluate the longitudinal length of the channel, 
the height profile, cross-sectioning the highest peak of the channel as well as the mica 
surface in its course, was used. 

Image processing. The averaged images of KcsA molecules were processed using a 
custom-made program running on Igor Pro 6.2. A reference molecule was picked up 
as a region of interest (ROI) and the molecules with similar features were collected 
using correlation pattern matching over the entire image. KcsA molecules with a 
normalized correlation coefficient over 0.9 to the reference image were averaged in 
four-fold symmetry. 

MD simulation. The MD simulations were performed with using the AMBER9 
package 41 . In order to minimize the dependency of the trajectories on the initial 
coordinates, ten different initial coordinates were generated for both the closed and 
open structures at the NVT ensemble (310 K). A periodic boundary condition was 
imposed. Electrostatic interactions were calculated using the particle mesh Ewald 
method 42 with a 0.8 nm cutoff in real space. The cutoff of the Lennard-Jones 
interaction was 0.8 nm. The bonds including those of H atoms, were constrained with 
the SHAKE algorithm 43 , enabling a time step of 2 fs. 

The system was prepared in the following manner. First, the dioleoylphosphati- 
dylcholine (DOPC) bilayer was equilibrated with the surrounding KC1 solution for 
10 ns. The empirical potentials were used for water (TIP3P) 44 , DOPC 45 , and ions 46 . 
The area per lipid after the equilibrium was confirmed to be the same as that reported. 
To avoid the aggregation of KC1 clusters 47 , parameters for ion developed by Dang et al. 
were used. The MD simulation was performed at the NPT ensemble (1 bar, 310 K) 



using Berendsen's barostat 48 . A vacant space was then created by removal of the 
DOPC and solution for filling the KcsA channel (PDB code: 1K4C) with F125 added 
to the C-terminal end (residues 22-125). The AMBER94 parameters were used for the 
KcsA channel 49 . E71 was protonated as in the previous study 50 . H25, El 18, El 20, and 
HI 24 were also protonated with the structure in the open state under an acidic 
condition (pH = 4.0). Equilibrium simulation was performed for 10 ns at NPT 
ensemble (1 bar, 3 10 K), where the C a in residues 22-124 were under the influence of 
a weak harmonic constraint (0.25 kcal/mol/A 2 ). The simulation for equilibration in 
the NVT ensemble (310 K) was continued for 10 ns after the constraint on the KcsA 
channel was removed. 

The coordinates of the C a in residues 30-114 were constrained to those of a 
1.7 nm-wide open structure (PDB code: 3F7Y) using a weak harmonic constraint 
(0.25 kcal/mol/A 2 ) for 50 ns. Next, the coordinates of the same C a were induced to 
those of a 2.3 nm-wide open structure (PDB code: 3F7V) with the same harmonic 
constraint for 50 ns. Finally, the coordinates of the same C a were constrained to those 
of 3.2 nm-wide open structure (PDB code: 3F5W) in the same way. Although the 
restraints were imposed on the C a atoms, the resultant structure was similar to that of 
3F7V in respect of the average C a - C a distances at T112 of the diagonal subunits 
(2.4 nm). 

For ten different initial coordinates, each coordinate was calculated for 20 ns under 
NVT ensemble without restraint and was used for the analyses. The total simulation 
time was 200 ns for both closed and open states. In the closed state, the system 
comprised of one KcsA channel, 119 DOPC, 30 K + , 58 CI", and 10704 water mole- 
cules. The system measured 77.3 X 80.5 X 87.1 A 3 . 

Single-channel current recordings. Planar lipid bilayer method was used for the 
single-channel current measurements as reported previously 38 . A Teflon sheet of 
0.5 mm thickness with a hole (<100 urn diameter) was used as a partition which 
divided the two (cis and trans) compartments of a chamber. The planar lipid bilayer 
was formed on a hole by painting 20 mg/mL rc-decane (Nacalai, Kyoto, Japan) 
solution of l-palmitoyl-2-oleoyl-s«-glycero-3-phospho-(r-rac-glycerol) (Avanti 
Polar Lipids, Alabaster, AL). Electrolyte solution in the cis and trans chambers 
contained 200 mM KC1 buffered by 10 mM HEPES (pH 7.5) and 10 mM succinic 
acid (pH 4.0), respectively. The reference electrode was placed in the cis 
compartment, to which KcsA channels reconstituted in proteoliposomes were added. 
The current was recorded through an Axopatch 200B amplifier and Digidata 1322 A 
digitizer (Molecular Devices) and the data were analyzed by using pCLAMP software 
(Molecular Devices, Sunnyvale, CA). The low pass filter was set to 2 kHz for the 
cut-off frequency, and the data were sampled at 5 kHz. 
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